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Non-extensive Boltzmann Equation and Hadronization 

T. S. Biro* and G. Purcsel^ 
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For a general, associative addition rule defining a non-extensive thermodynamics we construct 
the strict monotonic fuirction, which transforms it to a normal extensive quantity. We investigate 
the evolution of the one-particle distribution in the framework of a two-body Boltzmann equation 
supported with a non-extensive energy addition rule. An H-theorem can be proven for the extensive 
function of the non-extensive entropy. The equilibrium distribution is exponential in the extensive 
quantity associated to the original one-particle energy with the non-extensive addition rule. We 
propose that for describing the hadronization of quark matter non-extensive rules may apply. 
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There are several experimental evidences on power-law 
tailed statistical distributions of single particle energy, 
momenta or velocity. In particular hadron transverse 
momentum spectra at central rapidity, which stem from 
elementary particle and heavy ion collisionSjCan be well 
fitted by a formula reflecting TOT-scahng|^ Q, IE H; IM1^ 

119 : 

/(pt) ~ (1 + • (1) 

Interpreting these spectra as a distribution in the trans¬ 
verse directions at zero rapidity, the single particle energy 
is given by if = rriT = for a relativistic parti¬ 

cle with mass m (in units setting the lightspeed to c = 1). 
Amazingly this formula describes exactly the Tsallis dis¬ 
tribution f{E), which was obtained by using theoreti¬ 
cal arguments of thermodynamical nature0. Distribu¬ 
tions with power-law like tail, in particular the Tsallis 
distribution, can be seen in many areas where statisti¬ 
cal models apply^E 0; ^ ■ It 

was investigated as a generic feature in the framework 
of non-extensive thermodvua,mics [TE Efl . Tsallis has 
suggested an expression for the entropy, also encountered 
earlier by others[^l2^ , which would be an alternative to 
the Boltzmann formula. From this, with a canonical con¬ 
straint on the total energy, the distribution m can be de¬ 
rived. Without being able to exclude a non-equilibrium 
interpretation of the power-law tail, it is tempting to in¬ 
vestigate the possibility that some non-exponential spec¬ 
tra would be a result of a particular form of equilib¬ 
rium, featuring characteristics of a non-extensive ther¬ 
modynamics. 

In this paper we propose a possible way to understand 
power-law tailed energy distributions as equilibrium so¬ 
lutions to a slightly generalized two-body Boltzmann 
equation. More specifically we show that this two-body 
Boltzmann-equation allows for non-exponential station¬ 
ary single-particle distributions, if the two-body distribu¬ 


tion factorizes, but the two-body energy is not extensive. 
The Tsallis distribution is a special case thereof. The 
pair-energy, typical for generating this distribution, is a 
product of the single particle energies. 

It is a widespread belief that only the exponential dis¬ 
tribution can be the stationary solution to the Boltzmann 
equation, but this statement is true only with a few re¬ 
strictions: i) if the two-particle distributions factorize, 
ii) the two-particle energies are additive in the single¬ 
particle energies {E 12 = Ei -\- E 2 ) and iii) the collision 
rate is multilinear in the one-particle densities. 

A generalization of the original Boltzmann equation 
has been pioneered by Kaniadakisll^ considering a gen¬ 
eral, nonlinear density dependence of the collision rates. 
An "H” theorem for the particular Tsallis form of 
the collision rate has been derived by Lima, Silva and 
PlastinoliE. Here we follow another ansatz, we mod¬ 
ify the linear Boltzmann equation in the energy balance 
part only: Instead of requiring Ei -\- E 2 = E^ -\- E 4 in a 
I -k 2 ^ 3 -k 4 two-body collision we consider a general, 
not necessarily extensive, rule: 

h{E4,E2) = hiE3,E4). (2) 

It is physically sensible to choose the function h{x, y) 
symmetric and satisfying h{E, 0) = h{0, E) = E. 
Also, for applying the same rule for subsystems com¬ 
bined themselves of subsystems, associativity is required: 
h{h{x,y),z) = h{x,h{y,z)). This way the same rule ap¬ 
plies for the elementary two-particle system as for large 
subsystems in the thermodynamical limit. It is known 
that the general mathematical solution of the associativ¬ 
ity requirement is given by 

h{x,y) = X-'^ {X{x) + X{y)), (3) 

with Ar(0) = 0 and X{f) being a continuous, strict mono¬ 
tonic function!^. Composing the formula with the 
function X and taking the partial derivative with respect 
to j/ at j/ = 0 one obtains an ordinary differential equation 
for X{x) with the solution 


X{E) 


X'{Q) 


dx 

Jo f(x,0)- 
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Due to X{h{Ei, E 2 )) = X{Ei)+X{E 2 ), the quasi-energy 
X(E) is an additive quantity and the rule @ is equiva¬ 
lent to 


X{E^) + X{E2) = XiEs) + X{Ei). (5) 

Applying such a general energy addition rule the rate 
of change of the one-particle distribution is given by 

= [ W(P1>P2;P3,P4)A [f{p3)f{P4) - f{Pl)f{P2)]- 

^P2,P3,P4 

( 6 ) 

with the symmetric transition probability W and the con¬ 
straint 

A = 5'^{pi+p2-P3-P4) 5 {h{Ei,E2) - h{E3,Ei)). (7) 

In equilibrium the distributions depend on the phase 
space points through the energy variables only and the 
detailed balance principle requires 

f{E,)f{E2) = f{E3)f{Ei). ( 8 ) 

With the generalized constraint 0 this relation is satis¬ 
fied by 


f{E) = f{0)exp{-X{E)/T) (9) 

with l/T = -/'(0)//(0) and X{E) given by eq.0. In 
the extensive case h{x,y) = x + y leads to X{E) = E, 
for the Tsalis-type energy addition rule p^Eof . 

h{x,y) = X+ y + axy, (10) 

one obtains X{E) = i ln(l -|- aE) and 

/(A) = /(0)(l+aA)-'/“^. (11) 

Since the energy addition rule © conserves the quan¬ 
tity h{Ei,E 2 ) in a microcollision, the new energies af¬ 
ter the collision also lie on the h=constant line. Due to 
the additivity of the quasi-energy, X{E), the total sum 
Xtot = is a conserved quantity. This rule 

can be applied in numerical simulations, too. Since our 
present goal is to find the equilibrium only, we may as¬ 
sume constant transition probabilities. 

Fig n presents results of a simple test particle simula¬ 
tion with the rule (tTUll . We mostly started with a uni¬ 
form energy distribution between zero and Eq = \ with 
a fixed number of particles N = 10^ (red, full line) . 
The one-particle energy distribution evolves towards the 
well-known exponential curve for a = 0, shown in the 
left part of FigUl This snapshot was taken after 200 
two-body collisions per particle (blue, short dashed line). 
The analytical fit to this histogram is given by 2e“^/^ 
(with T = Eq/2 = 0.5 in this case). An intermediate 
stage of the evolution after 0.4 collisions per particle is 
also plotted in this figure (green, long dashed line). Using 
the prescription m with a = 1, the stationary solution 




FIG. 1: (Color online) Evolution towards the Boltzmann 
distribution for a = 0 (left part) and towards the Tsallis one 
for a = 1 (right part) using h(Ei,E 2 ) = Ei + E 2 + aEiE 2 . 
For detailed explanation see text. 


becomes a Tsallis distribution. The numerical evolution 
from the uniform energy distribution can be inspected 
in the right side of FigHI The fit to the final curve is 
given by f{E) = 2.6(1 -|- aE)~^'^. All distributions are 
normalized to one. 

It is in order to make some remark on the energy con¬ 
servation. For h{x,y) = a; -|- y we simulate a closed sys¬ 
tem with elastic collisions: The sum, U = F'i, does 
not change in any of the binary collisions. The situa¬ 
tion changes by using a non-extensive formula for h{x, y), 
like in the case of the Tsallis prescription. With a con¬ 
stant positive (negative) a, the bare energy sum, U, is 
decreasing (increasing) while approaching the stationary 
distribution. This is typical for open systems gaining or 
loosing energy during their evolution towards a station¬ 
ary state. 

One may incline to consider the conserved quasi¬ 
energy, X{E), as an in-medium one particle energy. 
The interesting point is that in general any prescrip¬ 
tion, h{Ei,E 2 ) - defining a version of the non-extensive 
thermodynamics is equivalent to considering a quasi¬ 
energy, X{E). For small energies one expects a restora¬ 
tion of the extensive rule and X{E) « E, A'(0) = 1. 
Whenever the pair energy is repulsive (attractive), 
h(^Ei, E 2 ) > El E 2 {h(Ei^ E 2 ) < El + E 2 ), a ris¬ 
ing quasi-energy is smaller (bigger) than the free one, 
x{E) < E {x{E) > f;)|33i. This leads to a tail of the 
stationary distribution in the free single particle energy, 
f{E), which is above (below) the exponential curve. This 
phenomenon is hard to distinguish from a power-law tail 
numerically!^ . 

The question arises that - constrained by the con¬ 
served number of particles, N = J fdT, and the total 
quasi-energy, Xtot = / f X{E) dT, - what is the proper 
formula for the entropy which grows when approaching 
the stationary distribution. If the addition rule of the 
non-extensive entropy, s, is given by hs{x,y), then the 
quasi-entropy, As(s), is additive and the total entropy is 
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Boltzmann entropy Tsallis entropy 




FIG. 2 : (Color online) Evolution of the Boltzmann and Tsal¬ 
lis entropies by applying the energy addition rule h{E\,E2) = 
El -\- E2 + aEiE2 with different values of the parameter a. 

given by 

XsiStot) = J fXsis{f))dr = Jaif)dr. (12) 

Its rate of change, Xs{Stot) = f f<^'{f)dr can be ex¬ 
pressed with the help of the Boltzmann equation I®. 

Assuming the symmetry properties 1 •*-> 2, 3 <--> 4 and 
(12) <-> (34) for the constrained rate factor ■u;i 234 = WA, 
one easily derives 

XsiStot) Stot = T / W;i234 ifsfi — / 1 / 2 ) X 

4 J 1234 

(^'(/i)+^'(/2)-a'(/3)-a'(/4)). (13) 

A definite sign for this quantity can be obtained due 
to the additivity of Xs{s{f)) = <j{fyf, leading to the 
unique solution Xsis{f)) = Blnf . For B = —ks 
(Boltzmann’s constant) Stot > 0 follows. Using the 
/cb = 1 unit system we arrive at Xs{sif)) = — In/, and 
the expression for the total additive quasi-entropy 
coincides with Boltzmann’s original suggestion. At the 
same time, applying a non-extensive addition rule for the 
entropy, hsix,y) = x + y + [1 — q)xy, as Tsallis did, we 
have As(s) = ln(l-|-(l—q)s) (Abe’s formula[2l|). and 

from Xsis) = — In / one obtains Tsallis’ entropy formula 

fi - f 

= irrf- ( 14 ) 

In numerical simulations we observe that the Tsallis-type 
non-extensive energy addition rule, applied to a test- 
particle simulation with two-body collisions, maximizes 
the Boltzmann-quasi-entropy (cf. Fig®. The station¬ 
ary distribution is nevertheless a Tsallis distribution (cf. 
Fig®. 

It is interesting to note that, as in many na.nersjT^. 
applying a canonical constraint on the bare total energy, 
U = JEfiE) dr, when seeking for the equilibrium distri¬ 
bution one obtains equivalent results for the Tsallis case. 


Instead of using the (due to the H-theorem guaranteed 
correct) formula, 

±ifXsisim = a + pXiE), (15) 

where Xs(s) and X{E) are both Tsallis-like expressions, 
the naive assumption of 

^ ifsif)) = a +PE, (16) 

also leads also to a power-law distribution. In the 
first case Xs{s{f)) = —In/ independently of the non- 
extensive rule definition for the entropy, and the equilib¬ 
rium solution is feq cx exp(—A(U)/T) = (1 -|- 
with P = 1/T. In the second case the particular form 
of s(/) is used and the non-extensive addition rule for 
the energy is ignored. (This case has nothing to do with 
our simulation.) The distribution resulting from lUTHl is 
feq CX (1 + aE/il + aaE))~^Sa'T) q = I — aT. 

In numerical simulations both the Boltzmann and the 
Tsallis entropy increase, disregarding the fluctuations 
due to the finite number of test particles, N = 10^, and 
due to finite binning resolution in energy, = 100 

while the maximal energy, Emax is changing). The in¬ 
crease is a trend, it is not rigorously fulfilled in each mi¬ 
croscopic collision; it is a common behavior of molecu¬ 
lar dynamical simulations. The ratio of the two entropy 
expressions also fluctuates somewhat, but the Tsallis en¬ 
tropy in the repulsive (attractive) case clearly stays big¬ 
ger (smaller) than the Boltzmann entropy. 

Physical realizations of non-extensive systems may be 
discovered depending on our knowledge about the micro¬ 
scopical forces influencing the particles during the pair- 
interactions. For such forces being repulsive, the canon¬ 
ical one-particle energy distribution has a tail above the 
exponential curve, for attractive interactions below. As a 
rule, as long as this modification is small, it is extremely 
difficult to see the non-exponential tail in the bare one- 
particle energy distribution both in experiments and in 
numerical simulations. Such power-law tails are promi¬ 
nent in elementary particle spectra in high energy experi¬ 
ments, but their traditional explanation does not assume 
an equilibrium state. 

In the quark-gluon plasma (QGP), or more generally 
in a parton matter before hadronization color non-singlet 
objects are the single particles. Eventually all form 
hadrons, in the soft sector perhaps by recombination and 
in the hard sector dominantly by fragmentation. In both 
cases a long-range interaction between color non-singlet 
partons, connected to the physical phenomenon confine¬ 
ment, is present in the background. In the following we 
consider a simple model for including this type of non- 
perturbative pair-interaction. 

For the sake of simplicity let us restrict ourselves to two 
body processes between color triplets and anti-triplets. 
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This is the most common way of meson formation. It 
is also an important part of baryon formation due to 
quark - diquark fusion. The pairs of such partons, while 
they constantly interact, are either in a color singlet or 
in a color octet state (in a QGP in one case from nine 
a singlet, otherwise an octet). The energy of the two- 
parton system is given by 

^ color state _ ^ E -h (l^) 

where the singlet channel should be attractive (relative 
to the free partons). The color average is supposed to 
be vanishing, = 0. This is certainly 

the case for interactions like in the Heisenberg-model 
of magnets, where the pair-potential is proportional to 
the product of symmetry generators in the correspond¬ 
ing spin representation. For SU(3) color this is also the 
case. The singlet charge is zero, the octet charge square 
is Qoctet = 3. The triplet and anti-triplet both have 
^triplet = 4/3. The Heisenberg-magnet-like interaction 
in color has therefore a factor of —8/3 for the singlet and 
a factor of 3 — 8/3 = 1/3 for the octet. Their degeneracy- 
weighted sum is zero. 

For considering the possibility of a non-Boltzmann dis¬ 
tribution in quark matter we further assume a Coulomb¬ 
like interaction. In this case from the 

binding in the color singlet channel. For the search af¬ 
ter a stationary single-quark distribution of the two-body 
Boltzmann equation in the octet channel it accounts to 
consider. 


qoctet ^ El + E2 + (18) 

The rest is kinematical consideration. We assume the 
coalescence of two massless partons to a (nearly) massless 
hadron. Due to the triangle inequality, the kinetic energy 
of the relative motion of two massless partons is non¬ 
negative. 


F; -l-kin. = + p- I - \P1 1 - |P2 I > 0. (19) 


For small relative angles between the momentum vectors, 
■i?, this is approximated by 


E 


rel.kin. 

12 


2E1E2 sin 2 (t?/ 2 ) 

El E2 


( 20 ) 


The sum of the individual parton energies in the same 
approximation is close to Ei + E 2 « P = |pi -I- P 2 1 • Even 
very hard partons with a high value of the total pair 
momentum, P, need a little relative motion for interact¬ 
ing: in the singlet channel to eventually form hadrons, 
in the octet channel to maintain a single-particle quark- 
distribution typical for the pre-hadronic phase. The sta¬ 
tionary version of this distribution, while detailed bal¬ 
ance is satished on the two-body level, is often found to 
be close to the Tsallis distribution. We propose that the 


above mechanism, from the comparison of eqs. (Cnii, di) 
and iinj leading to 


Pc = 1/a 


2P 

sm^{d/2y 


( 21 ) 


may be in the background of such findings. Asymptoti¬ 
cal freedom is recovered as for very fast partons Ec —>^ 00 
with P —> 00 , and so the one particle energies of a col¬ 
liding pair become additive. 

In conclusion we have investigated deterministic, non- 
extensive energy addition rules in two-body collisions. 
We have pointed out that instead of the one-particle en¬ 
ergy a quasi-energy is conserved by such rules in each 
collision, leading to a non-Boltzmannian stationary dis¬ 
tribution in the bare one-particle energy. In particular 
the Tsallis distribution is obtained by using a Tsallis- 
type non-extensive energy addition rule. The corre¬ 
sponding conserv ed q uasi-energy is identical to that pro¬ 
posed by Q. Wangle . Modifications to the extensive en¬ 
ergy addition rule may have to be considered if there is 
a statistically important pair-interaction between parti¬ 
cles. The stationary canonical distribution becomes ex¬ 
ponential in the conserved quasi-energy X{E), but it is 
non-exponential in terms of the free particle energy E. 
The Boltzmann entropy, Sb = Xs{Stot) = —f /In/dF, 
is never decreasing and reaches its maximum at this 
distribution. Alternative expressions for the entropy, 
in particular the one promoted by Tsallis, correspond 
to a non-extensive entropy addition rule which defines 
Xs(s). Notably, in the Tsallis case also the naive in¬ 
tegral, St = f /s(/)dr, seem to increase in numerical 
simulations as the elementary collisions proceed. The 
Tsallis-type energy and entropy addition rules are corre¬ 
lated by g = 1 — aT. 

As a possible physical realization we have proposed a 
mechanism leading to nearly Tsallis-distributed quarks in 
quark matter and hadrons which eventually form. This 
mechanism considers a color state dependent pair-energy 
based on general arguments describing a quantum sym¬ 
metry. The essential clue leading to our result then hides 
in the use of a virial theorem which connects the color 
interaction with kinematical factors of the quark pair. In 
a certain approximation the modification of the familiar 
two-body energy conservation factor in the Boltzmann 
equation receives a term proportional to the product of 
single-quark kinetic energies to leading order in the ul- 
trarelativistic expansion. The Tsallis distribution turns 
out to be an approximation next simplest to the original 
Boltzmannian oue|33||. 

Our picture naturally connects the processes maintain¬ 
ing a possible stationary distribution among colored par- 
tons with the hadronization process. In the approxima¬ 
tion discussed in this paper the power of the power-law 
tail in hadronic spectra equals to the power occurring 
in the single-quark distribution in quark matter. As a 
consequence mesonic and baryonic powers are also equal 
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to each other. This agrees with experimental findings 
well, although the recombination assumption predicting 
a baryonic to mesonic power ratio of 3 : 2 also cannot be 
excluded, when considering the relatively high error bars. 
This result is, however, simpler and more generic. Here 
only a balance between kinetic and potential energy in 
the relative motion of quarks has been assumed besides 
some basic color properties of the pairwise interaction. 
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